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£j ' ABSTRACT 

Context. It is widely accepted that the large obliquity of Uranus is the result of a great tangential collision (GC) with an Earth size 
proto-planet at the end of the accretion process. The impulse imparted by the GC had affected the Uranian satellite system. Very 
recently, nine irregular satellites (irregulars) have been discovered around Uranus. Their orbital and physical properties, in particular 
those of the irregular Prospero, set constraints on the GC scenario. 

Aims. We attempt to set constraints on the GC scenario as the cause of Uranus' obliquity as well as on the mechanisms able to give 
origin to the Uranian irregulars. 

Methods. Different capture mechanisms for irregulars operate at different stages on the giant planets formation process. The mech- 
anisms able to capture the uranian irregulars before and after the GC are analysed. Assuming that they were captured before the GC, 
we calculate the orbital transfer of the nine irregulars by the impulse imparted by the GC. If their orbital transfer results dynamically 
implausible, they should have originated after the GC. We then investigate and discuss the dissipative mechanisms able to operate 
later. 

Results. Very few transfers exist for five of the irregulars, which makes their existence before the GC hardly expected. In particular 
Prospero could not exist at the time of the GC. Different capture mechanisms for Prospero after the GC are investigated. Gas drag 
by Uranus 'envelope and pull-down capture are not plausible mechanisms. Capture of Prospero through a collisionless interaction 
seems to be difficult. The GC itself provides a mechanism of permanent capture. However, the capture of Prospero by the GC is a low 
probable event. Catastrophic collisions could be a possible mechanism for the birth of Prospero and the other irregulars after the GC. 
Orbital and physical clusterings should then be expected. 

Conclusions. Either Prospero had to originate after the GC or the GC did not occur. In the former case, the mechanism for the 
origin of Prospero after the GC remains an open question. An observing program able to look for dynamical and physical families is 
mandatory. In the latter case, another theory to account for Uranus' obliquity and the formation of the Uranian regular satellites on 
the equatorial plane of the planet would be needed. 

Key words. Planets and satellites: general - Planets and satellites: formation- Solar System: general- Solar System: formation 



1. Introduction 

Very recently, rich systems of irregular satellites (hereafter irregulars) of the giant planets have been discovered. Enabled by the 
use of large-format digital images on ground-based telescopes, new observational data have increased the known population of 
Jovian irregulars to 55 (Sheppard et al. 2003), the Saturnian population to 38 (Gladman et al. 2001, Sheppard et al. 2005a 2006a) 
and the Neptunian population to 7 (Holman et al. 2004 Sheppard et al. 2006b). The Uranian system is of particular interest since a 
population of 9 irregulars (named Caliban, Sycorax, Prospero, Setebos, Stephano, Trinculo, S/2001U2: XXIV Ferdinand, S/2001U3: 
XXII Francisco and S/2003U 3: XXIII Margaret) has been discovered around Uranus (Gladman et al. [19981 120001 Kavelaars et al. 
2004 Sheppard et al. 2005b). The discovery of these objects provides a unique window on processes operating in the young Solar 
System. In the particular case of Uranus, their existence may cast light on the mechanism responsible for its peculiar rotation axis 
(Parisi & Brunini [19971 Brunini et al.|20Tj2l (here after BP02)). 

Send offprint requests to: M. Gabriela Parisi 
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Irregulars of giant planets are characterized by eccentric, highly tilted with respect of the parent planet equatorial plane, and in 
some case retrograde, orbits. These objects cannot have formed by circumplanetary accretion as the regular satellites but they are 
likely products of an early capture of primordial objects from heliocentric orbits, probably in association with planet formation itself 
(Jewitt & Sheppard 2005 ). It is possible for an object circling about the sun to be temporarily trapped by a planet. In terms of the 
classical three-body problem this type of capture can occur when the object passes through the interior Lagrangian point, L2, with a 
very low relative velocity. But, without any other mechanism, such a capture is not permanent and the objects will eventually return 
to a solar orbit after several or several hundred orbital periods. To turn a temporary capture into a permanent one requires a source 
of orbital energy dissipation and that particles could remain inside the Hill sphere long enough for the capture to be effective. 

Although currently giant planets have no efficient mechanism of energy dissipation for permanent capture, at their formation 
epoch several mechanisms may have operated: 1) gas drag in the solar nebula or in an extended, primordial planetary atmosphere or 
in a circumplanetary disk (Pollack et al l 19791 Cuk & Burns I20031 . 2) pull-down capture caused by the mass growth and/or orbital 
expansion of the planet which expands its Hill sphere (Brunini |l995l Heppenheimer & Porco [T977] >, 3) collisionless interactions 
between a massive planetary satellite and guest bodies (Tsui 1999 ) or between the planet and a binary object (Agnor & Hamilton 
2006), and 4) collisional interaction between two planetesimals passing near the planet or between a planetesimal and a regular 
satellite. This last mechanism, the so called break-up process, leads to the formation of dynamical groupings (e.g. Colombo & 
Franklin 1971, Nesvorny et al. [2004). After a break-up the resulting fragments of each progenitor would form a population of 
irregulars with similar surface composition, i.e. similar colors, and irregular shapes, i.e. light-curves of wide amplitude. Significant 
fluctuations in the light-curves of Caliban (Maris et al. 2001) and Prospero (Maris et al. I2007at and the time dependence observed 
in the spectrum of Sycorax (Romon et al. 120011 ) suggest the idea of a break-up process for the origin of these bodies. 

Several theories to account for the large obliquity of Uranus have been proposed. Kubo-Oka & Nakazawa ( 11995ft , investigated 
the tidal evolution of satellite orbits and examined the possibility that the orbital decay of a retrograde satellite leads to the large 
obliquity of Uranus, but the large mass required for the hypothetical satellite makes this possibility very implausible. An asymmetric 
infall or torques from nearby mass concentrations during the collapse of the molecular cloud core leading to the formation of the 
Solar System, could twist the total angular momentum vector of the planetary system. This twist could generate the obliquities of 
the outer planets (Tremaine |1991| ). This model has the disadvantages that the outer planets must form before the infall is complete 
and that the conditions for the event that would produce the twist are rather strict. The model itself is difficult to be quantitatively 
tested. Tsiganis et al. (|2p05 ) proposed that the current orbital architecture of the outer Solar System could have been produced from 
an initially compact configuration with Jupiter and Saturn crossing the 2:1 orbital resonance by divergent migration. The crossing 
led to close encounters among the giant planets, producing large orbital eccentricities and inclinations which were subsequently 
damped to the current value by gravitational interactions with planetesimals. The obliquity changes due to the change in the orbital 
inclinations. Since the inclinations are damped by planetesimals interactions on timescales much shorter than the timescales for 
precession due to the torques from the Sun, especially for Uranus and Neptune, the obliquity returns to small values if it is small 
before the encounters (Hoi et al. 120071 ). 

Large stochastic impacts at the last stage of the planetary formation process have been proposed as the possible cause of the 
planetary obliquities (e.g. Safronov |1969l ). The large obliquity of Uranus (98°) is usually attributed to a great tangential collision 
(GC) between the planet and an Earth-size planetesimal occurred at the end of the epoch of accretion (e.g., Parisi & Brunini 1997 
Korycansky et al.Q990). The collision imparts an impulse to Uranus and allows preexisting satellites of the planet to change their 
orbits. Irregulars on orbits with too large semimajor axis escape from the system (Parisi & Brunini 1997), while irregulars with 
a smaller semimajor axis may be pushed to outer or inner orbits acquiring greater or lower eccentricities depending on the initial 
orbital elements, the geometry of the impact and the satellite position at the moment of impact. The orbits excited by this perturbation 
must be consistent with the present orbital configuration of the Uranian irregulars (BP02). 

In an attempt to clarify the origin of Uranus obliquity and of its irregulars, we are using in this study the most updated information 
on their orbital and physical properties. 

In Section 2, we improve the model developed in BP02 for the five Uranian irregulars known at that epoch and extend our study 
to the new four Uranian irregulars recently discovered by Kavelaars et al. (2004) and Sheppard et al. (2005b). The origin of these 
objects after the GC is discussed in Section 3, where several mechanisms for the origin of Prospero are investigated. The discussion 
of the results and the conclusions are presented in Section 4. 



2. Transfer of the irregulars to their current orbits: 

Assuming the GC scenario, the transfers of the nine known irregulars to their current orbits are computed following the procedure 
developed in BP02 for the five irregulars known in 2002. We present improved calculations using a more realistic code to compute 
the evolution of the irregulars current orbital eccentricities. 

If the large obliquity of Uranus has been the result of a giant tangential impact, the orbits of preexisting satellites changed due 
to the impulse imparted to the planet by the collision. The angular momentum and impulse transfer to the Uranian system at impact 
were modeled using the Uranus present day rotational and orbital properties as imput parameters (BP02). 

Just before the GC, the square of the orbital velocity vi of a preexisting satellite of negligible mass is given by: 




(1) 



M. Gabriela Parisi et al.: Constraints to Uranus' Great Collision IV 



3 



r being the position of the satellite on its orbit at the moment of the GC, a\ its orbital semiaxis and my the mass of Uranus before 
the impact. The impactor mass is m, and G is the gravitational constant. After the GC, the satellite is transferred to another orbit 
with semiaxis «2 acquiring the following square of the velocity: 

2 / 2 1 \ 

v 2 = Gimu + w,) . (2) 

\r a 2 ) 

We set v\ — A v 2 and v\ — B (1 + mt/rrtu) v^, where A and B are arbitrary coefficients (0 < A < 1, B > 0), v e being the escape 
velocity at r before the GC. 

The semiaxis of the satellite orbit before (a\) and after (ai) the GC verify the following simple relations: 

r r 

fli = , Gb = • (3) 

2(1 -A) 2(1 -E) 

If A < B then a\ < In the special case of B — 1, the orbits are unbound from the system. If A > B then a\ > the initial 
orbit is transferred to an inner orbit. When A = B, the orbital semiaxis remains unchanged {a\ = ai). 

The position r of the satellite on its orbit at the epoch of the impact may be expressed in the following form: 



2 G my 
r " (AV)2 



B'-A 



VAcos^i y/(B' - A) + A cos 2¥ 



2, (4) 



with B' — B (1 + nii/mu). Since stochastic processes can only take place at very late stages in the history of planetary accretion (e.g. 
Lissauer & Safronov |1991| ), the GC is assumed to occur at the end of Uranus formation (e.g. Korycansky et al. 1990). The mass of 
Uranus after the GC, (m, .+mu), is taken as Uranus' present mass. *P is the angle between v\ and the orbital velocity change imparted 
to Uranus AV. An analytical expression for AV is derived in BP02 assuming that the impact is inelastic (Korycansky et al. 1 19901 ) as 
a function of m,, the impact parameter of the collision b, the present rotation angular velocity of Uranus Q, the spin angular velocity 
which Uranus would have today if the collision had not occurred SIq. and a which is the angle between and 



2Rl 
AV = — - 

5b 



2QQn cos a 

£r + 



(l + a.)2 ( i + "L)4 (i + + j^)2 

v my 7 v imu ytijj 7 v im\j ' 



1/2 



(5) 



Ru being the present equatorial radius of Uranus. A collision with the core itself was necessary in order to impart the required 
additional mass and angular momentum (Korycansky et al. 119901 . Since b is an unknown quantity, we take its most probable value: 
b={2/3)Rc, where Rc is the core radius of Uranus at the moment of collision assumed to be 1.8 x 10 4 km (Korycansky et al. 1990 
Bodenheimer & Pollack 1986). The results have a smooth dependence with the impactor mass which allows us to take m, ~ lm ffl 
(Parisi & Brunini [l997] l. 

The minimum eccentricity of the orbits before the collision is given by: 

e lmi „ = 2 (1 - A) - 1 if A < 0.5 , e lmi „ = 1 - 2 (1 - A) if A > 0.5, (6) 
while the minimum eccentricity of the orbits after the collision is: 

eimin = 2 (1 - fi) - 1 if B<0.5 , e 2mi „ = 1-2(1 - B) if B > 0.5. (7) 

The minimum possible value of AV (AV m ,„) is obtained from Eq. © for an initial period To- 20 hrs (Tq = 2n/Qo) and a = 70° 
(BP02). Therefore, although simulations of solid accretion produce in general random spin orientations (e.g. Chambers 2001), we 
further assume To- 20 hrs and a-lQ° in order to set a maximum bound on Eq. ©. Upper bounds in a\ [a\u) and a 2 («2m) are 
obtained from Eq. (O through Eqs. and (0 taking AV=AV„„„ with V F=180°, i.e., assuming the impact in the direction opposite 
to the orbital motion of the satellites and taking the positive sign of the square root in Eq. (HJ: 

Gmu(B'-A)2 _ G m v (B' - A)2 

d\M - ^= pr- , a 2 M - 1= p— ■ ■ \o) 

(A V,„,„)2 ( 1 - A)( VB 7 - VA) 2 (A V mi „)2 ( 1 - B)( Vfi 7 - VA) 2 

For each A, we calculate the value of B (B= B' /(l + nti/mu)) corresponding to the transfer to aiM- a, where a is the present 
orbital semiaxis of each one of the Uranian irregulars shown in the second column of Table 1 in units of Ru- From Eq. (0, this 
value of B provides the minimum possible value of e2„,i„, e-i m , that the orbit of each irregular may acquire at impact for each initial 
condition A and every initial condition for To, a, *P and m,, i.e., if a transfer of a given orbit (A,B) is not possible for ¥=180°, To= 
20 hrs and a = 70°, the same transfer (A,B) is not possible for any other incident direction of the impactor and for any other value 
of 7q and a either. 

Since the orbits of the irregulars are time dependent, the orbital evolution of the five Uranian irregulars known in 2002 was 
computed in PB02 by numerical integration of the equations of the elliptical restricted three body problem formed by the Sun, 
Uranus and the satellite. In this paper, we present the orbital evolution of the nine known Uranian irregulars for 10 5 yrs using the 
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Satellite 


r s [km] 


a[Ru] 


ineaii 


ai[Ru] 


e, 


(Aa)/a,) 


(Ae/e,) 




Caliban 


49 


283 


0.191 


287.5 


0.1973 


1.602 x 10- 2 


3.289 x 10 


-2 


Sycorax 


95 


482 


0.541 


485 


0.5436 


6.307 x 10~ J 


4.795 x 10 


-3 


Prospero 


15 


645 


0.432 


648 


0.4342 


4.585 x 10- j 


4.997 x 10 


-3 


Setebos 


15 


694 


0.581 


701.8 


0.5853 


1.123 x 10~ 2 


7.366 x 10 


-3 


Stephano 


10 


314 


0.251 


333 


0.2781 


6.058 x 10~- 


0.1080 




Trinculo 


5 


336 


0.218 


361.6 


0.2501 


7.623 x 10~- 


0.1475 




Ferdinand 


6 


813 


0.660 


839.3 


0.6701 


3.241 x 10~ 2 


1.533x 10" 


2 


Francisco 


6 


169 


0.142 


280.6 


0.3593 


0.6604 


1.530 




Margaret 


5.5 


579 


0.633 


649 


0.6705 


0.1209 


5.917 x 10 


-1 



Table 1. Present parameters of the Uranian irregulars and orbital damping due to gas drag exerted by Uranus extended envelope. 
r s and a are the present physical radius and the present orbital semiaxis of the irregulars. e mean is their calculated mean eccentricity 
tabulated in Table 2. a, and e, are the orbital semiaxis and eccentricity just after the GC, while (Aa)/a,) and (Ae/e,) are the damping 
of these orbital elements since the epoch of the GC until the contraction of Uranus envelope. 



Satellite 


&mean 


^max 


^min 


Caliban 


0.191 


0.315 


0.072 


Sycorax 


0.514 


0.594 


0.438 


Prospero 


0.432 


0.571 


0.305 


Setebos 


0.581 


0.704 


0.463 


Stephano 


0.251 


0.381 


0.121 


Trinculo 


0.218 


0.237 


0.200 


Ferdinand 


0.660 


0.970 


0.393 


Francisco 


0.142 


0.187 


0.093 


Margaret 


0.633 


0.854 


0.430 



Table 2. Variation of the eccentricity of the Uranian irregulars due to Solar and giant planet perturbations over a period of 10 5 yrs. 



symplectic integrator of Wisdom & Holman ( 119911 1, where the perturbations of the Sun, Jupiter, Saturn and Neptune are included. 
The mean (e mean ), maximum (e max ) and minimum (e„„„) eccentricities are shown in Table 2 for all the known Uranian irregulars. 

The transfer of each satellite from its original orbit to the present one is possible only for those values of (A,B) which satisfy 
the condition e2„, < e„ wx . A satellite did not exist before the impact if it has no transfer and the satellites with the widest range of 
transfers are those with the highest probability of existing before the impact. 

The transfers within a range of 20 Ru around each present satellite semiaxis («2M- fl ± 20 Ru', a taken from Table 1) for all the 
Uranian irregulars are shown in Fig. 1 . There are few transfers for Setebos, Ferdinand and Margaret. This makes the existence of 
these satellites before collision little probable. The only transfers for Trinculo and Prospero are close to the pericenter of an eccentric 
initial outer orbit (e [m > 0.58 for Prospero and e lm > 0.62 for Trinculo). The minimum eccentricity after collision e-i m for Trinculo is 
in the range [0. 16-0.23], very close to e max (0.237). This result gives a very low probability for the existence of Trinculo before the 
GC. For Prospero e2 m is in the range [0.52-0.57], e2 m ~ e max (0.571). Therefore this satellite could not exist before the GC. If the 
present large obliquity of Uranus was caused by a large impact at the end of its formation, Prospero had to originate after the event. 
Relating the origin of the outer Uranian system to a common formation process, all the Uranian irregulars probably were originated 
after the GC. The possible post-GC origin of Prospero and the other Uranian irregulars is discussed in the following section. 

3. Origin of Prospero after the Great Collision 

In this section, we analyze the possibility that Prospero have been captured after the GC. We investigate the possible dissipative 
mechanisms able to produce its permanent capture taking into account that the giant impact is assumed to have occurred at late 
stages in the planetary accretion process. 

3.1. Gas drag by Uranus' envelope 

Bodenheimer & Pollack (1986) and Pollack et al. {1996) studied the formation of the giant planets by accretion of solids and gas. 
In their model, the so called core instability scenario, when the mass of the core of the planet has grown enough a gaseous envelope 
begins to form around it. For Uranus, its envelope extended until its accretion radius which was ~ 500 Ru at the end of Uranus' 
formation (Bodenheimer & Pollack 1986). The formation of Uranus is completed when there is no more nebular gas to accrete. 
Otherwise, gas accretion by proto-Uranus would have continued towards the runaway gas accretion phase and the planet would 
have now a massive gaseous envelope. Bodenheimer & Pollack (1986) obtained that after the end of accretion, the radius of the 
envelope of proto-Uranus remained almost constant (~ 500 Ru) over a time scale of 10 4 yrs and then contracted rapidly to ~ 8 Ru 
in 10 5 yrs. The final contraction to the present-day planetary radius occurred on a slower timescale of 10 8 yrs. 

Korycansky et al. (1990) carried out hydrodynamical calculations of the GC for a large set of initial conditions at the end 
of accretion. They found a sharp transition between the cases where almost all the mass of the envelope of Uranus remained 
after the impact and those where it was almost entirely dispersed by the impact. This implies that the impact should have not 
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Fig. 1. The transfers capable of producing the present orbits of the Uranian Irregulars. A (B) is the square of the ratio of the satellite's 
speed just before (after) the impact to the escape velocity at the satellite's location just before (after) the impact. e\ m (e2 m ) is the 
minimum eccentricity of the orbits before (after) collision. The full-black line A=B divides the upper region (the current orbits 
arise from inner orbits) and the lower region (the current orbits arise from outer orbits). The value of e max tabulated in Table 2 is 
shown on the dashed line for comparation with e^. Trinculo: empty down triangles, Caliban: full circles, Sycorax: empty rhombus, 
Ferdinand: empty squares, Francisco: full down triangles, Margaret: empty up triangles, Prospero: full up triangles, Setebos: empty 
hexagons, Stephano: empty circles. 



dispersed the envelope, as there would have been no nebular gas to re-accrete on the planet. They showed that the envelope reacts 
hydrodynamically at impact and it expands outward. After the shock the gas falls back on the core over a timescale of few hours 
being the final result a readjustment instead of a catastrophic transformation. The timescale for this hydrodynamical process is much 
shorter than the orbital period of the irregulars which is of the order of years. We may then assume that the GC did not change the 
envelope density profile. 
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The extended envelope of Uranus could be in principle, a source of gas allowing the capture of Prospero and the other irregulars 
after the GC. Assuming that the GC did not change the envelope profile, we fit from Fig.l of Korycansky et al. ( 1990) the density 
profile of Uranus' gaseous envelope before the GC, p g = c te R~ 4 g cm -3 with c te = 10 36 and R being measured in cm. It gives a 
nebular density of ~ 4 x 1CT 13 g cirT 3 at the boundary of 500 Ru in agreement with the minimum mass nebula model. 

As a first approximation, we compute the ratio of gas mass traversed by a body of density p s and radius r s in a characteristic 
orbital period P, to the mass of the body. Assuming for the body a circular orbit of radius R, we calculate the so-called /3 parameter 
(Pollack et al. \WM : 

P 2nRp„nr 2 3np„R 

P=~= 4 3 ~ W 

t ±np s ri 2p s r s 

where t is the characteristic timescale for changing any of the orbital parameters. For the permanent capture to occur /3 cannot be 
very small, /3 > 0.04 (Pollack et al. H979t . Using Eq.© and assuming a nebula density 10 times that of the minimum nebula model 
(c te = 10 37 ) for an object the size of Prospero (p s = 1 .5 g cm 4 ) and at Prospero's pericenter (R= 278 Ru), /3 ~ 9 x 10~ 5 , which is too 
small to affect the orbit of Prospero. 

Following BP02, we now investigate with more detail the possible effect of gas drag on the Uranian irregulars after the GC due 
to Uranus' extended envelope before its contraction to its present state. Following the procedure of Adachi et al. (119761 1, we obtain 
the time variations of the eccentricity e and semiaxis a of each Uranian irregular. The drag force per unit mass is expressed in the 
form: 



where v re i is the relative velocity of the satellite with respect to the gas. In computing the satellite mass m, a satellite mean density 
p s of 1.5 g cirT 3 is taken for all the satellites (http://ssd.jpl.nasa.gov/7sat_phys_par). The drag coefficient Co is ~ 1 and r s is each 
satellite radius taken for each satellite from Table 1. 

Assuming that the orbital elements are constant within one Keplerian period (the variations of a and e are very small), we 
consider the rates of change of the elements averaged over one period, that is: 



Since after the end of accretion the gas density in the outer regions of the envelope contracts rapidly, we have integrated Eqs.([PTTi 
back in time on 10 4 yrs for the 9 Uranian irregulars. We have taken v re / as the satellite orbital velocity since we have assumed a 
null gas velocity. This assumption maximizes the orbital damping for retrograde satellites which allows us to set upper bounds in 
the damping effect for the orbital eccentricity and semiaxis of all the retrograde irregulars. The mean eccentricity e mean and the 
actual semiaxis a from Table 1 were taken as the initial conditions for the integrations. We take p g - c te R 4 g cirT 3 with c te = 
10 36 and R = a(l - e 2 )/(l + ecos#), a being measured in cm. The orbital damping is shown in Table 1, where a,- and e, are the 
initial semiaxis and eccentricity just after the GC at the end of accretion and Aa=(a,- - a) and A<?= (<?,■ - e mean ), are the damping in 
the orbital semiaxis and eccentricity. Stephano, Trinculo and Margaret had experienced little orbital evolution while Francisco had 
suffered a large orbital damping (see Table 1). The permitted transfers of Fig. 1 would increase for Francisco since the condition eim 
< e, should be then satisfied. However, the orbital evolution of the larger satellites, in particular of Prospero, results negligible. Even 
increasing the nebula density by a factor of 10 (c te =10 37 ) the damping of the orbital elements of Prospero is too small with Aa/a,= 
0.046 and Ae/e,= 0.048. This satellite had not experienced any orbital evolution due to Uranus' extended envelope and then could 
not have been capture by gas drag after the GC. 

The pressure forces acting on a body traveling through the gas not only decelerates it, but also subjects it to stresses. If the 
stress is greater than the strength of the body, the body is fractured in fragments of different size. The fragments move away one 
another since drag forces vary inversely with size and act to separate them. The average pressure on the forward hemisphere of a 
non-rotating, spherical body as it moves through the gas with relative velocity v re i is approximately equal to the dynamic pressure, 
Pdyn = \Pg v2 re i (Pollack et al j 19791 >. The body will fragment into pieces if pd yn > Q, where Q is the compressive strength. Values 
of Q on the order of 3 x 10 6 dyne crrr 2 are needed to shatter strong (e.g. rock/ice) targets ( which is 10 times lower than the value 
adopted for asteroids), while compressive strength on the order of 3 x 10 4 dyne cirT 2 are appropriate for relatively weak (snow-like) 
targets (Farinella & Davis IT9961 Stern[1996). For a body on circular orbit at the present pericenter of Prospero (R= 278 Ru), where 
Vrei is the circular speed around Uranus at R and taking c te = 10 37 for p g , pd yn ~ 0. 17 dyne cirT 2 , and at Sycorax's pericenter pd yn ~ 1 
dyne cm 2 . In both cases pd yn « Q and Prospero could not have been originated by the dynamical rupture of a parent object. 

A collision may fracture the parent body but if the energy at impact is not sufficient to disperse the fragments, drag forces may 
act to separate them against their mutual attraction. The relative importance of these effects is measured by the ratio ej of the drag 
force on a given fragment j to the gravitational force acting on j by the other fragments i (Pollack et a! J19791 >: 



F = -Cp g v 2 reI ,C = 



C D nr 2 s 



(10) 



2m 




(ID 




(12) 
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where F ci j is the gravitational force between the particles i and j, and F D j is the drag force on the particle j of radius rf 

C D p g nr 2 .v 2 

F Dj = — '—, (13) 

The fragment j is dispersed by the gas if ey > 1 (Pollack et al !19791 l. In order to estimate the order of magnitude of the effect, we 
computed Eq.dTSb for j=Prospero at Prospero's pericenter taking into account the gravitational attraction due to another fragment of 
equal size, and at Sycorax pericenter taking into account the gravitational force of Sycorax on Prospero. In the first case, we obtain 
€j= 10 and in the second case at Sycorax pericenter ej= 0.04. We then conclude that pressure forces are not strong enough neither 
to fragment a possible parent object from which Prospero originated, nor to disperse its fragments. In addition Prospero suffered no 
orbital evolution due to gas drag and could not have been capture by Uranus 'envelope after the GC. 



3.2. Pull-down capture 

Within the GC scenario, runaway of the cores of the planets occurred during the first stages of accretion but stopped for each 
embryo after it reached a size of about 1000 Km. At 10-35 AU the final mass distribution contained several hundreds of Mars-size 
(or larger) bodies dominating the mass of the residual disk. Beauge et al. (2002), investigated the effects of the post-formation 
planetary migration on satellites orbits. They obtained that if the large-body component (composed of Mars-size bodies) dominated 
the mass of the residual disk, the presently accepted change in the orbit of Uranus of ~ 3 AU is too large and it is not compatible 
with the observed distribution of its satellites. Even an orbital change of ~ 1.5 AU already causes sufficient instabilities to eject all 
the Uranian irregulars. Pull-down capture caused by the orbital expansion of the planet could then not be a plausible mechanism for 
the origin of Prospero and the other irregulars. Pull-down capture caused by the mass growth of the planet after the GC would not 
be possible given the impact is assumed to have occurred at the end of the accretion process when there was no more mass to be 
accreted by the planet. 



3.3. Collisionless interactions 

Within the framework of the restricted three-body problem, a capture is always followed by an escape. To end up with a long term 
capture, the satellite has to dissipate energy in a short time. The entrance energy AE within the gravitational field of the planet is 
(Tsui [1999l l: 

A£ = -2.15// 2/3 (l-5)^, (14) 
a p 

yi- M P /M Q and 5 « 1, where M p and a p are the mass and orbital semiaxis of the planet . 

Tsui (I1999K suggested a permanent capture mechanism where a guest satellite encounters some existing inner orbit massive 
planetary satellite causing its velocity vector to be deflected keeping the irregular in orbit around the planet. In this way, the 
effective two-body potential would be about twice the entrance energy AE of the guest satellite. The radius Ri of the orbit of the 
guest satellite after deflection is then given by: 

0.23 

Ri = —/ a v (15) 

In the case of Uranus, for a minimum entrance energy of 5-0, the minimum permanent orbital radius of the guest satellite is 
Ri = 955 Ry - This value of Ri is much larger than the present semiaxis of Prospero (see Table 1), making the capture of Prospero by 
this mechanism implausible. 

The fact that binaries have recently been discovered in nearly all the solar system's small-body reservoirs suggests that binary- 
planet gravitational encounters could bring a possible mechanism for irregulars capture (Agnor & Hamilton 2006). One possible 
outcome of gravitational encounters between a binary system and a planet is an exchange reaction, where one member of the binary 
is expelled and the other remains bound to the planet. Tsui ( 1 999 ) extended the scenario of large angle satellite-satellite scattering to 
the formation of the Pluto-Charon pair assuming that Pluto was a satellite of Neptune and that Charon was a guest satellite. Through 
Eqs. dT4l > and ( fT5b . the conditions for the escape of the pair was found. Following their scenario, let us consider the hypothesis that 
Prospero was a member of a guest binary entering Uranus' field, with energy density A£),,„, above the minimum density given by 
Eq.dT4b and 5=0. A close encounter with Uranus could result in disruption of the binary, leading to the ejection of one member 
and capture of the other. The minimum semiaxis R\ is given by Eq.dl3b. However, even this scenario seems to be unlikely since the 
semiaxis of Prospero is smaller than 955 Ru- 



3.4. Collisional interactions: Break-up processes 

Collisional interactions between two planetesimals passing near the planet or between a planetesimal and a regular satellite, the so 
called break-up process, leads to the formation of dynamical groupings (e.g. Colombo & Franklin 1971, Nesvorny et al. 2004). The 
resulting fragments of each progenitor body after a break-up will form a population of irregulars expected to have similar surface 
composition, i.e. similar colors, and irregular shapes, i.e. large temporal variations in the light curve as these irregular bodies rotate. 
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The critical rotation period (T c ) at which centripetal acceleration equals gravitational acceleration for a rotating spherical object 

is: 

< 3n ^ 



Tc = \ 7 ^\ , (16) 

where G is the gravitational constant and p„b is the density of the object. With p o /,=0.5, 1, 1.5 and 2 g cirT 3 , T c = 4.7, 3.3, 2.7 and 2.3 
hrs, respectively. The rotation period of Prospero is about 4 hrs (Maris et al. 2007a), and it seems unlikely that light and dark surface 
markings on a spherical Prospero could be responsible for its light curve amplitude of 0.2 mag (Maris et al. l2007ab . Even at longer 
periods, real bodies will suffer centripetal deformation into aspherical shapes. For a given density and specific angular momentum 
(H), the nature of the deformation depends on the strength of the object. In the limiting case of a strengthless (fluid) body, the 
equilibrium shapes have been well studied (Chandrasekhar [l987l ). For H < 0.304 [in units of (GM 3 /?^) 1 ' 2 ], where M (kg) is the 
mass of the object and R sp he is the radius of an equal-volume sphere, the equilibrium shapes are the oblate Maclaurin spheroids. 
For 0.304 < H < 0.390 the equilibrium figures are triaxial Jacobi ellipsoids. Strengthless objects with H > 0.390 are rotationally 
unstable to fission. The Kuiper Belt objects (KBOs), being composed of solid matter, clearly cannot be strengthless. However, it is 
likely that the interior structures of these bodies have been repeatedly fractured by impacts, and that their mechanical response to 
applied rotational stress is approximately fluid-like. Such "rubble pile" structure has long been studied in the asteroid belt (Farinella 
et al fTMTIl a nd in the Kuiper Disk (Sheppard & Jewitt 120021 Jewitt & Sheppard [20021 Romanishin & Tegler [19991 ) . Farinella & 
Davis d 19961 ) obtained that KBOs larger than about 100 km in diameter are massive enough to survive collisional disruption over 
the age of the solar system, but may nevertheless have been internally fractured into rubble piles. 

Whether a collision between an impactor and a target results in growth or erosion depends primarily on the energy of the impact 
and the mass and strength of the target. If the mass of the impactor is small compared to the mass of the target m s , the energy 
required at impact to result in a break-up is given by: 

1 , 3 Gm] 

- msVcol > ms S + -—, (17) 

where v co ; is the collision speed, S is the impact strength, R ms is the radius of the target and y is a parameter which specifies the 
fraction of collisional kinetic energy that goes into fragment kinetic energy and is estimated to be ~ 0.1 (Farinella & Davis 1996). 
The speed of the fragments is critical when the target has a gravity field. Fragments moving slower than the local escape speed 
re-accumulate to form rubber pile structures. 

We attempt to investigate whether Prospero could be a collisional fragment or if a collision on a primary Prospero would result 
in a rubber pile structure. 

In computing Eq.dT7t, v 2 col = v 2 e +v 2 in ^ where v e is the scape speed at the target surface and v,„y is the typical approach velocity 
of the two objects at a distance large compared with the Hill sphere of the target. For two bodies colliding in the Kuiper disk, v,„/ is 
given by (Lissauer & Stewart [19931 ): 

^/^ 2 (r 2+ 4 (18) 

where is the keplerian velocity, e is the mean orbital eccentricity and i the mean orbital inclination of the KBOs. We take (e)—2(i) 
( Stern [19961 ) . Eq.dTTb was computed using Eq.([T8b for values of e in the range [0.01-0.1] and orbital semiaxes of the KBOs in the 
range [30-60] AU. In computing the mass of the target m s , we consider the radius R ms of the KBOs in the range [10- 500] km and 
densities between [0.5 - 2] g cirT 3 . Impact strengths in the range [3 xl0 4 -3 xlO 6 ] erg cm ~ 3 were taken. We obtain that targets with 
Rms < 210 km suffer disruption for all the values of these parameters in the Kuiper Disk. Prospero has a radius of just only 15 km. If 
Prospero originated from the Kuiper Belt, it would have been more likely a collisional fragment rather than primary body. Prospero 
would preserve an irregular shape after disruption since it is such small object that is unable to turn spherical because its gravity 
cannot overcome material strength. Prospero could have been captured during the break-up event if the two KBOs collided within 
Uranus' Hill sphere, which could be possible for a minimum orbital eccentricity of the original KBOs of 0.37. 

We also consider the case in which the target is a satellite of Uranus which collides with a KBO that enters the Hill sphere of 
the planet. Eq.(fT7b remains valid but the following expression of v,„/ is considered: 



Vinf 



(19) 



where v lp is the velocity of the KBO with respect to Uranus and v sp is the satellite orbital velocity, at the epoch of th event. We assume 
that the satellite orbital velocity is circular. In order to get bounds in the relative velocity, we take two values of v;„/, v- ln f-v\ p + v sp . 
For Vi P , we assume that the KBO describes a hyperbolic orbit around Uranus during the approach giving: 

2 2G(mi+mu) GM 

V inf = + ' ( 20 ) 

a s dkb 

where we assumed GM^/a^b as the relative velocity between the KBO and Uranus far from the encounter, (m, + ma) is the present 
mass of Uranus and a s the orbital semiaxis of the satellite. We calculate Eq.dTTl) using Eqs.(fT9ll and ( f20b for in the range [20,60] 
AU and a s [100 -700] Ru for the same values of S, R ms and densities we have used for the collisions among KBOs. We obtain that 
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for all the possible parameters, any Uranus' satellite with radius R ms < 1000 km suffers disruption if it collides with a KBO larger 
than 10 km. This process would lead to the formation of two clusters of irregulars, one associated to the preexisting satellite and 
the other to the primary KBO. This process has the disadvantage that it is unlikely that the preexisting satellite were formed from a 
circumplanetary disk as regular satellites given the large orbital semiaxis required for this object. 

Break-up processes predict orbital clustering. However, no obvious dynamical groupings are observed at the irregulars of Uranus. 
A further intensive search of more faint irregulars around Uranus is needed in order to look for dynamical and physical families. 

3.5. The CG itself as a possible capture mechanism 

We now turn to the question of whether the GC itself could have provided a capture mechanism (BP02). Since all the transfers with 
A > B' lead to a more bound orbit, this process might transform a temporary capture into a permanent one (see Section 2 and Fig. 1). 
Moreover, a permanent capture could even occur from an heliocentric orbit (transfers with A-l). 

It is interesting to estimate the number of objects N in heliocentric orbits at the time of the GC, at distances from Uranus less 
than or equal to 300 Ru- Assume that the GC occurred when Uranus was almost fully formed, meaning that its feeding zone was 
already depleted of primordial planetesimals. We assume that the objects passing near Uranus at that time were mainly escapees 
from the Kuiper belt. Using the impact rate onto Uranus and the distributions of velocities and diameters given by Levison et al. 
(2000), and assuming that the mass in the transneptunian region at the end of the Solar System formation was 10 times its present 
mass, a back-of-the-envelope calculation gives one object of diameter D > 20 km passing at a distance R < 300 Ru from Uranus 
every 6 yrs at the end of accretion (BP02). The typical crossing time Tc among protoplanets in the outer Solar System is larger than 
one millon years (Zhou et al. 2007 ). The number of objects passing near Uranus during a timescale Tc is then 167000, which gives 
a probability of 6 x 10~ 6 for the capture of an object at about 300 Ru by the GC. This low rate of incoming objects, makes the 
possibility of the capture of all the irregulars from heliocentric orbits difficult. Even the capture of a single object, Prospero (note 
that Prospero could not have an orbit bound to the planet before the GC), turns out to be low probable. Since temporary capture can 
lengthen the time which a passing body can spend near the planet, a more plausible situation arises if we assume that the GC could 
produce the permanent capture of one or more parent objects which were orbiting temporarily around Uranus being the present 
irregulars the result of a collisional break-up occurring after the GC. 



4. Discussion and Conclusions 

It is usually believed that the large obliquity of Uranus is the result of a great tangential collision (GC) with an Earth-sized proto- 
planet at the end of the accretion process. We have calculated the transfer of angular momentum and impulse at impact and have 
shown that the GC had strongly affected the orbits of Uranian satellites. We calculate the transfer of the orbits of the nine known 
Uranian irregulars by the GC. Very few transfers exist for five of the nine irregulars, making their existence before the GC hardly 
expected. In particular, Prospero could not exist at the time of the GC. Then, either Prospero had to originate after the GC or the 
GC did not occur, in which case another theory able to explain Uranus' obliquity and the formation of the Uranian regular satellites 
would be needed. It is usually believed that the regular satellites of Uranus have accreted from material placed into orbit by the GC 
(Stevenson et al. |1986b . 

Within the GC scenario, several possible mechanisms for the capture of Prospero after the GC were investigated. If the Uranian 
irregulars belong to individual captures and relating the origin of the outer uranian system to a common formation process, gas 
drag by Uranus' envelope and pull-down capture seem to be implausible. Three-body gravitational encounters might be a source 
of permanent capture. However, we found that the minimum permanent orbital radius of a guest satellite of Uranus is ~ 955 Ru 
while the current semiaxis of Prospero is 645 Ru- The GC itself could provide a mechanism of permanent capture and the capture 
of Prospero could have occurred from a heliocentric orbit as is required within the GC scenario, but due to the low rate of incoming 
objects it turns out to be difficult. Break-up processes could be the mechanism for the origin of Prospero and the other irregulars 
in the frame of different scenarios. Prospero might be a fragment of a primary KBO fractured by a collision with another KBO. 
The fragment could have been captured by Uranus if the two KBOs had a minimum orbital eccentricity of 0.37. Prospero could be 
a secondary member of a collisional family originated by the collision between another satellite of Uranus and a KBO where the 
parent satellite of Prospero could have been captured by any mechanism before or after the GC. This process has the disadvantage 
that it is unlikely that the preexisting satellite were formed from a circumplanetary disk as regular satellites given the large orbital 
semiaxis required for this object. Since collisional scenarios require in general high collision rates, perhaps the irregulars were 
originally much more numerous than now. Then, Prospero and also the other irregulars might be the result of mutual collisions 
among hypothetical preexisting irregulars (Nesvorny et al. 2003 2007) which could have been captured by any other mechanism 
before the GC. 

The knowledge of the size and shape distribution of irregulars is important to know their relation to the precursor Kuiper Belt 
population. It could bring valuable clues to investigate if they are collisional fragments from break-up processes occuring at the 
Kuiper Belt and thus has nothing to do with how they were individually captured later by the planet, or if they are collisional 
fragments produced during or after the capture event (Nesvorny et al. 2003, 2007). The differential size distribution of the Uranian 
irregulars approximates a power law with an exponent g=1.8 (Sheppard et al. 2005b). If we assume that the size distribution of the 
nine irregulars with radii greater than 7 km extends down to radii of about 1 km, we would expect about 75 irregulars of this size or 
larger (Sheppard et al 20050). 

The nuclei of Jupiter family comets are widely considered to be kilometer- sized fragments produced collisionally in the Kuiper 
Belt (Farinella & Davis 1996). Jewitt et al. (2003 ) compared the shape distribution of cometary nuclei in the Jupiter family with the 
shape distribution of small main-belt asteroids of similar size (1 km -10 km) and with the shape distribution of fragments produced 
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in laboratory impact experiments. They found that while the asteroids and laboratory impact fragments show similar distribution 
of axis ratio ((b/a) ~ 0.7), cometary nuclei are more elongated ((b/a) ~ 0.6). They predict that if comets reflect their collisional 
origin in the Kuiper Belt followed by sublimation-driven mass loss once inside the orbit of Jupiter, small KBOs should have average 
shapes consistent with those of collisionally produced fragments (i.e., (b/a) ~ 0.7). To date, constraints on the shapes of only the 
largest KBOs are available. Prospero being a bit larger than cometary nuclei, displays a variability of 0.2 1 mag in the R band (Maris 
et al. [2007 a |. This corresponds to an axis ratio projected into the plane of the sky, b/a of 0.8. The knowledge of the size and shape 
distribution of irregulars would shed light in the size and shape distribution of small KBOs as well as on the irregulars capture 
mechanism. 

Colors are an important diagnostic tool in attempting to unveil the physical status and the origin of the Uranian irregulars. 
In particular it would be interesting to assess whether it is possible to define subclasses of irregulars just looking at colors, and 
comparing colors of these bodies with colors of minor bodies in the outer Solar System. Avalilable literature data show a dispersion 
in the published values larger than quoted errors for each Uranian irregular (Maris et al. [2007a, and references therein). We have 
concluded in Maris et al. (2007a), that the Uranian irregulars are slightly red but they are not as red as the reddest KBOs. 

An intensive search for fainter irregulars and a long term program of observations able to recover in a self consistent manner 
light-curves, colors and phase effects informations is mandatory. 
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